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^ i 1 Introduction 

^^ 

*^t ' The double heavy mesons, i.e. the heavy quarkonia {ccj, (bb) and Be meson 

as well as its excited states, are non-relativistic bound-state. Their static 

properties and most decays may be computed quite precisely by means of the 

^^ , potential model (PM) inspired in QCD and relevant theories as well, whereas 

0^ ' their hadronic production was not understood so well, although the hadronic 

(— I , production of J/tp is so important, not only for perturbative QCD (PQCD) 

O ' theory itself but also for its applications, for instance, it is used as a probing 

tool to detect specific signature in nuclear processes such as that of quark 

gluon plasma (QGP) etc. Substantial progress on the hadronic production has 

(~| , been achieved in recent years. Last year, the discovery of the double heavy 

meson Be opened up new challenges for studying the problems. Moreover, at 

tree level the Be meson, being a ground state and in flavor non-singlet, decays 

by weak interaction only, and the branching ratios of the decays are sizable^ 

Hence in addition to the mesons B and D, we will be able to study the weak 

decays of the heavy flavors 6-quark and c-quark with the meson Be- 

2 Advances and Problems 

Fresh progressions on the double heavy mesons are outlined below. 
2.1 Advance on Fragmentation Functions of the Double Heavy Mesons 



"A talk presented in the OCPA Conference "Recent Advances and Cross-Century Outlooks 
in Physics" , March 18-20, 1999, at Atlanta 

''On the contrary, the ground states of heavy quarkonia (in flavor singlet) decay through 
strong and electromagnetic interactions mainly, that most of the branching ratios of the 
possible weak decays are so tiny beyond the capacity of the present experimental techniques, 
although the widths of the weak decays have the same order as those of the meson Be- 
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Figure 1: Factorization of an Inclusive Production for Hadron Collisions and Definition of 
Fragmentation Function (F.F.). 



It is known that^ perturbative QCD (PQCD) works efficiently with its 
factorization theorenu. Namely the 'hard part' involved a concerned process 
can be calculated perturbatively, but the 'soft part', e.g. structure functions 
(S.F.s) and fragmentation functions (F.F.s), being in non-perturbative nature, 
relating to the 'color confinement' and independent on the specific process, are 
just treated as 'universal factors'. The factorization theorem and definition of 
a fragmentation function may be illustrated precisely by Fig.P and (El) with an 
inclusive process. 

da = "^ dxi dx2 / dx^F^^ixi, hf) •i^4(^2,MF) 

■ daij^k,x{xi,X2.,xz,^iF) ■ Dl{xy,,^F), (1) 



where F^^ (xi, /if) and F^ {x2, IJ-f) are the S.F.s to depict the probabilities for 
finding a parton i in the hadron hi with momentum fraction xi and a parton j 
in the hadron /12 with momentum fraction X2 respectively; daij_^k,x (2^1, X2, x^, ^f) 
is the cross-section corresponding to the 'hard subprocess' (the SP-part in 
Fig.|l|); the factor D^{x^,^f) is just the F.F., which depicts a probability for 
the produced parton k to form a hadron h with a momentum fraction 2:3 of the 
parton k. Here /if, denoting an energy scale, is noted where all the factors in 
the equation are matched, and sometimes it is omited for simplification if there 
is no cofusion. The fragmentation functions, being considered independent of 
the specific process, contain nonperturbative effects, that generally they are 
considered not to be calculable, but may be determined by measurement (s). 
Thus F.F.s are important blocks for PQCD calculations. 

Recently, an advance is to be realized that the fragmentation functions 
(F.F.s) for double heavy mesons can be factorized further, and some of them. 
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Figure 2: The Further Factorization for a Fragmentation Function of the Double Heavy 
Meson Be (an example). 



for the so-called colaij|Singlet ones, even can be calculated out precisely without 
any free parameteD0a. Taking a F.F. of the meson Be as an example, let me 
with Fig.g to illustrate the 'further factorization'. Note here that the virtual 
gluon which creates a pair of c, c as depicted in Fig. 2 is crucial for the further 
factorization. It is because that the gluon carries a time-like momentum, and 
always has q^ > Arn^ ^ ^qcd ' ^^ ^^^^ ^^ i^ guaranteed the factor involving the 
gluon being perturbative. Namely, due to the very virtual gluon, the further 
factorization makes sense, and we may calculate out the perturbative part of 
the F.F. precisely. In the cases, when the 'residual' factor of the further factor- 
ization is in color singlet, although being non-perturbative, the residual factor 
may be directly related to the wave functions of the potential model directly, 
i.e. it can be calculated precisely without any free parameter. Therefore a 
fragmentation function for a double heavy meson can be calculated out either 
totally if the residual factor is in color-singlet which may be related to the wave 
function of potential model directly, or partly i.e. only the virtual gluon part 
may be calculated out and the residual factor, as a factor to be determined, is 
left in the 'final result' if the residual factor is in color non-singlet. 

To take an example and to see the general features, a typical fragmentation 
function, calculated out by the method firsfl, has the behavior: 



D^^iz,mBj^ai{m')\iP{0)\ 



^z{l^zf 



{[2aiz — 3(a2 — ai)(l - a2z) 



{a2Z-lf 
(2 - z)](l - a2z)z + 6(1 + aizf{l - a2z) - 801022^(1 - z)}, (2) 



where |'0(O)| is wave function at original of the Be meson; a'^{m'^ ) is the QCD 
couphng constant at m^ ; and ai — mc/mB^, 0,2 — mb/mB^- 

2.2 Advances: 'Puzzle' and Solution for Charmonium Production 

Thanks to the reahzation about the further factorization for the frag- 
mentation functions of double heavy mesons, various fragmentation functions 
not only for the mesons Be, B*, ■ ■ ■ but also for the heavy quarkonia (cc) and 
(bb) have been computed. For the hadronic charmonium production, on one 
hand, the experimental advance is the data have been accumulated for many 
years, and new achievements in experimental techniques have made it possible 
to separate the prompt production events of charmomia J/?/;, V'',Xc, • • • from 
those events as decay-products of the produced B mesons in hadron collisions 
successfully; on the other hand, the theoretical advance on the fragmentation 
function computation for those of color singlet: the fragmentation functions 
have been calculated out precisely by means of the wave functions obtained by 
potential model (all the parameters appearing in the calculation having been 
determined well elsewhere) , therefore thorough comparisons between theoreti- 
cal predictions and experimental data become possible and interesting. Indeed 
the comparisons were carried through just after the advances being achieved. 

An important result of the comparison, the so-called 'V'' surplus puzzle', 
rose and confirmed. Namely, of the prompt production of J/ip and ip' in hadron 
collisions, the experimental data are greater than theoretical prediction by one 
to two orders in magnitude and the behavior of the Pt (the transverse momen- 
tum of the produced mesons) dependence for experimental data and theoretical 
prediction is also inconsistent, if only the color singlet F.F. is taken into ac- 
count. We would like to note here that the advance on the computation of 
the fragmentation functions for double heavy mesons is crucial to recognize 
the puzzle, namely, without the advance, the fragmentation functions would 
be able only to be determined by experimental measurement in whole, that 
the deviation (the puzzle) would be 'absorbed' into the determination to dis- 
appear. Based on analysis of the puzzle, Braaten and Fleming proposed the 
so-called color octet mechanism firscl. Following them, many groups apnhed 
the proposition to various processes widely, including to nuclear processe£l to 
explore the consequences of the new mechanism. Since in the proposition a new 
undetermined parameter, the so-called color-octet matrix element (C.O.M.E.) 
which is of non-perturbative feature and uncalculable, is introduced, hence the 
solution of the 'surplus puzzle' essentially is to determine the matrix element 
first and then to test it with the other available data. Only with constant ma- 
trix elements determined by fitting data to obtain correct behavior of the Pt 
dependence of the production has 'absolute meaning' for the proposition. The 
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Figure 3: 1. The Color Octet Component cc{[^5i,8]). 



referencoS paid attention to the Pt dependence of the production on large Pt 
mainly, but afterwards others tried to have a complete determination of the 
matrix elements for various color octet components by fitting whole bdjia.vior 
of PT-dependence (at large and small Pt) of the hadronic productioiErullJ. 

Taking the charmonium meson J/'ijj as an example let me outline the key 
point of the proposition. A physical state of J/ip, in general, should be ex- 
pressed in Fock space: 



I J/^) - O(z;0)|cc(p5il])) + Oiv)\cci['Pj8])g) + Oiv')\cc{['So8])g) 
+ 0{v^)\cc{['Si8])gg) + 0{v^)\cci['Djl, 8])gg) + • • ■ , 



(3) 



where v is 4- velocity of the constituents of the meson J/V', and, due to J/ip 
being heavy nonrelativistic bound state, v^ ~ 0.3 is acceptable; g indicates a 
'valance' gluon; and the factor 0{v^) at the front of each term indicates the 
relative order of the term which is obtained based i-on the counting rules of 
the effective theory: nonrelativistic QCD (NRQCD^^. It is known that J/^p 
is a 'flavor singlet' meson, thus it has the same quantum number in flavor as 
a gluon, especially, its component in the Fock space expansion (ccpS'iS]) has 
the same quantum numbers of a gluon. As a result of the fact, J/?/; production 
in hadron collisions will present a quite complicated feature i.e. in a sense it 
may be considered as a kind of 'mixing' between the components of higher 
Fock space with a gluon. The fragmentation function of a gluon to J/ 4', 
Dg'^{z,fj,F = 2mc), may be used as an example to illustrate how to solve the 
surplus by the 'mixing'. For charmonium, in general, we have v'^ ~ 0.3 and 
as(4TO^) ~ 0.2 ~ 0.3, therefore, there are components for the fragmentation 
function, which have contribution in the same order of magnitude, according 
to naive NRQCD power-counting rule, are listed below: 

1. The color octet component (ccpS'i,8]) (Figure 3:1): 



D'^/^{z, p5i,8]) = {Oi^''['S^,8]) X r,^(,,ps„8])(^). 



(4) 
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Figure 3: 2. The Color Octet Component cc(pPj,8]). 
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Figure 3: 3. The Color Singlet Component cc([^Si,l]). 



2. The color octet component (ccpP,/,8]) and {cc[^So,S]) (Figure 3:2): 

D'J^iz, ['Pj,8]) = {Oi/'^['Pj,8]) X r,^(,,-[3p,,8]),(z), (5) 

D'J^iz, [i5o,8]) = (Oi/^['S„,8]) X r,_,(,,[.So,8]),(^). (6) 

Note that for shortening in Fig. 3:2, only {cc['^Pj, 8]) is presented, but {cc[^So, 8]) 
being similar, is not. 

3. The color singlet component (ccp5i, 1]) (Figure 3:3): 



'>J/->I'[3 



Di'nz, [^5i, 1]) = H^n^^i, 1]) X r,^cc-(ps.a]),,(^). 



(7) 



Based on power counting and as indicated in Figs. 3:1-3 no matter the com- 
ponent in the color-octet or in color-singlet, each may contribute to the F.F. 
in the same order of magnitude, that we should treat them 'equally'. 

2.3 The Problem,(s) of the Proposition 

Indeed the color octet proposition is very attracting because it is on NRQCD 
formalism and can solve the puzzle by choosing the color octet matrix elements 
(C.O.M.E.) properly. Whereas the C.O.M.E. cannot be calculated due to its 
nonperturbative nature (there is a long way to go also for Lattice QCD), and 
furthermore, exactly to say, the proposition has not been well-proved yet. 
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The proof of the proposition should be either to show the determination of 
C.O.M.E. is independent of the chosen processes, or one and more of its char- 
acteristic signatures are observed in experiments. Whereas the status of this 
kind of proof now is that a lot of suggestions on the characteristic signatures 
are proposed but no one has been observed in experiments definitely|j; for the 
determination of C.O.M.E., the situation is neitipr definite: An example for 
the situation of the determination of C.O.M.Effl'EB is shown in Table 1. 

Table 1: Color-octet matrix element linear combinations. 



Matrix elements 


(0^''*(35i))(mGey3) 


(O-V,^, , <°«'y^»»(.nGey3) 


NRQCD scaliiiH order 

Cho et ain 

Beneke et Epi 

Sridhar et alt3 


(6.6 ±2.1) X 10-3 
(1.06 ±0.14) X 10-2 
(1.26 ±0.33) X 10-2 


(2.2 ±0.5) X 10-2 
(3.14 ±0.58) X 10-2 



Note: those of Beneke et atl are taken only for (fcy) = 0. The values in the 
table are obtained only based on the data of Tevatron, and those from the 
lepton(photon)-production and fixed target ones are not involved because of 
more complications. The values of C.O.M.E. in the table are quite 'scattering' 
that I think it is hard to draw any conclusion about proving the proposition. 

2.4 Be Mesons: Theoretical Predictions and Experimental Discovery 

The meson Be is the unique one of the possible double heavy mesons 
in fiavor non-singlet. Being a double heavy meson, its properties may be 
estimated well with the potential model and recent advances. 

Based on the potential model—the estimated mass of the meson is not very 
large in the region 6.2 — 6.4 Ge\0. 

The lifetime of the meson may be estimated by the spectator mechanism 
and the 'CKM favor' annihilation of the two components c-quark and 5-quark 
into c, s-quarks and r^ , j/-leptons (there is helicity suppression that those of 
annihilating into the other possible quark pairs and lepton pairs besides into 
c, s-quarks and t+, j/-leptons may be ignored safely). The two spectator terms 
are the same as those for B meson and for D meson respectively, thus they 
may relate to the widths of the B meson and the D meson i.e. we may estimate 
the contribution of the two spectator terms to the width of the meson Be from 
the experimental measurements on the lifetimes of the mesons B and D'^{Ds) 
respectively, and the annihilation contribution may be computed easily with 

■^One clean signature of the proppaition was proposed in referencit3, but considering to the 
backgrounds pointed out in ReftJ and experimental errors, one cannot conclude a clean 
signature has been observed. The other suggestions are worse than the mentioned one. 



the wave function (decay constant) o£_thij|_meson Be- The esthxiated hfetime 
of the meson B^ is about 0.4 x lO^i^yyy 

Since the mass of the meson is not very great, the production of the meson 
had not been expected so 'difficult' as reahzed now: Of the present facihties, 
for experimental detecting only at Tevatron, enough number of Be mesons can 
be produced; in e~^e~ colliders only at Z resonance LEP-I marginal number of 
the mesons can be produced; of the planned facilities, at such a high 
hadronic collider LHC, numerous mesons will be able to be producec 
Quantitative values depend on the cuts of transverse momentum Pt and ra- 
pidity y quite sensitively. Comparative study on the production calculated 
by full a^ calculation of PQCD and by simplified 'fragmentation mechanism' 
showed that due to convolution in the total crosa,section the differences of the 
two approaches were smoothed out accidentalljcS. 

As for the weak decays of the meson Be, there are many decay channels 
that contain a charmonium in the final state. Especially, when the charmo- 
nium in final state is J/^, the decays may be used as a characteristic signature 
for identifying the produced meson Be because of the detecting advantages of 
J/ip' narrow width and sizable branchingja. tin t p a lepton pair II etc. To com- 
pute all the decays, there are many wayEa't-lE3, but let nie-ipick up one, the 
so-called generalized instantaneous approximation approaclu3i22l to illustrate a 
little more in detail. The non-relativistic wave functions of the mesons Be and 
J/ip may be calculated out in potential model framework, but there may be a 
great (relativistic) momentum recoil (ms^ — 6.3 GeV, mji^ ~ 3.0 GeV, i.e. 
the produced J/iji may move relativistically in C.M.S. of Be). The so-called 
generalized instantaneous approximation approach is to start with a relativistic 
formulation, which is based on Bathe-Salpeter (B.S.) equation for the double 
heavy mesons and Mandelstam method for the transition matrix elementc3, 
to write down the relevant decay matrix element and then to make an 'instan- 
taneous approximation' on the whole matrix element in a similar way as the 
'original instantaneous approximation' on B.S. equation done by Salpeter first. 
Finally as the result of the approach, the matrix element is related to certain 
operators sandwiched by the nonrelativistic wave functions (the instantaneous 
'limit' of the B.S. equations) properly. Let me quote the semileptonic decay 
widths estimated by the approach and some other approaches in Table 2. 

Table 2: Exclusive semileptonic decay width for various modes 



A 



RefercnckJ (in lO"" eV) 



Mode 



Be -+ ric + 1+VI 
Bc^J/tp + l+ui 



Rcfcrciicda (in 10"" eV) 



14.2 
34.4 



10.6 
38.5 



The meson Be was searched for at LEP-I, but not seeiB. It was just 

8 



discovered by CDF Collaboration last yeaE3. Now no one doubts the discovery, 
although CDF Collaboration was based on about 20 events of the cascade 
decays Be — > J/ip + ^^i J/i^ ~^ /^"""M" ^-nd reported three measurements with 
quite large error bars only. The three are its mass, lifetime and a combination 
of the production cross section and the decay branching ratio of the meson. 
In summery for comparison, the CDF results of the three measurements: 

ms^ =6.40 ±0.39 ±0.13, (8) 

TB^ = OAGt^ollistat.) ± Om{syst.) x lO^^^^^ (g) 

and a combination of the production cross-section and decay branching ratio 
'normalized' by those of B meson's: 

a- BrjB, -^ J/^jj + 1+ + ■ ■ .) 

aBr{B+ -^ J/^ + K+) 

^ 0.132+°;°^^(stai.) ± 0.031(syst.)+°:o2o('*/e**"^e). (10) 

The ratio i?, the paass tob^ and the lifetime tb^ all are in the region of theo- 
retical prediction^. 

Since the rest decay channels have not been observed yet, thus to shorten 
the paper, I would not list them here, but just to mention the most interesting 
one Be — > J/^tt: roughly speaking, the width is small: just about one tenth 
of Be -^ J/ip + l^i^i, so we may expect it will be observed soon in RUN-II 
of Tevatron. Since each particle of the final state in the decay can be well- 
detected, so it is the best channel to measure the mass of the meson Be- 

3 Outlooks 

The advances and problem(s) on double heavy mesons have been outlined. 
It certainly is breakthrough progress, that of a physical state not only the color- 
singlet components but also the color-octet components in Fock space Eq.(3) 
may play very substantial roles, although it still needs to be further proved. 
The progress should have wide consequences to be explored. Reviewing the 
advances on the topics, many subjects are just at beginning i.e. there are many 
'things' needed to be study further. 

RUN-II of Tevatron will present its the first results within two or three 
years. In the next run: RUN-II, the total statistics will be raised by one or 
more order of magnitude, and the detectors will also be improved greatly so the 
systematic errors will be suppressed and the ratio of the signal to background 
will be raised much etc. In addition, considering to LHC to be built, we believe 
that in the very near future on the color-octet mechanisms of heavy quarko- 
nium production, systematic studies of the meson Be and much deeper under- 
standing must be achieved. Namely motivated by the foreseeing experimental 



achievements, including searching for the signatures of the color-octet proposi- 
tion, the C.O.M.E. measurements, accumulation of more Be events and more 
accurate theoretical calculations on Be meson etc must gain great progress 
soon. 

One specific point which I would like to address here is that a comparative 
study of the hadronic production of the heavy quarkonia and the meson Be 
is very helpful to clarify up the situation of the color-octet proposition. If 
the proposition is a correct solution for the original 'surplus' puzzle then it 
is necessary that for the production of Be meson there is no the 'surplus' 
at all, because the meson Be is in flavor non-singlet that the color-singlet 
component contribution in the production is always to play a dominant role, 
on contrary, a heavy quarkonium in flavor singlet, as indicated by Figs. 3:1-3, 
can gain the 'surplus' by the 'mixing' with a gluon i.e. a color-octet component 
can be produced in a lower order of PQCD than the color-singlet one although 
the former has a smaller possibility to form the hadron finally (C.O.M.E. is 
smaller). 
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